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Abstract 

The Lifshitz-type formulas describing the free energy and the force of the van der Waals interac- 
tion between an atom (molecule) and a single-wall carbon nanotube are obtained. The single-wall 
nanotube is considered as a cylindrical sheet carrying a two-dimensional free electron gas with 
appropriate boundary conditions on the electromagnetic field. The obtained formulas are used 
to calculate the van der Waals free energy and force between a hydrogen atom (molecule) and 
single-wall carbon nanotubes of different radia. Comparison studies of the van der Waals inter- 
action of hydrogen atoms with single- and multi-wall carbon nanotubes show that depending on 
atom-nanotube separation distance the idealization of graphite dielectric permittivity is already 
applicable to nanotubes with only two or three walls. 

PACS numbers: 73.22.-f, 34.50.Dy, 12.20.Ds 
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I. INTRODUCTION 



Van der Waals forces play a dominant role in the interaction of atoms and molecules with 
carbon nanostructures at separations larger than about one nanometer. These forces are of 
fluctuation nature and are described in second order perturbation theory as a dipole-dipole 
interaction. Theoretical description of the van der Waals forces as a quantum phenomenon 
was pioneered by London. The general theory of the van der Waals forces in the framework of 
quantum statistical physics was developed by Lifshitz jl, 2] using the concept of a frequency- 
dependent dielectric permittivity. Lifshitz's theory describes both the interaction between 
two semispaces or plates of finite thickness separated by a gap and between a microparticle 
(an atom or a molecule) and a flat surface of a macrobody. In the latter case it is commonly 
referred to as the Casimir- Polder force {sj]. 

An understanding of the underlying mechanisms of microparticle — nanostructure inter- 
action is of importance for many applications including the problem of hydrogen storage 
in carbon nanostructures 4|. However, the application of the Lifshitz theory to the lat- 
ter case encounters serious problems because the boundary surface of carbon nanotubes is 
not flat and single-wall nanotubes are not characterized by the macroscopic concept of a 
frequency-dependent dielectric permittivity. Because of this, most theoretical work on the 
van der Waals interaction in layered structures and between hydrogen atoms (molecules) 
and a graphite sheet or carbon nanotubes was done using the phenomenological density 



functional theory (see, e.g., Refs 



10 



llj). Some other approximate approaches 



were also used, e.g., the nonrelativistic perturbation theory for degenerate levels of a two- 



level atomic system 12], ^M- ^ * s known, however, that in some cases the approximate 



and phenomeno 



interaction 



0, 



ogical approaches do not provide a precise description of the van der Waals 



151 ]. This is true [16|, for instance, when one uses the density functional 
theory (especially with linear-density approximation). 

Recently, the scope of the Lifshitz theory of the van der Waals forces was much widened 
by successful application to the interpretation of precision measurements of the Casimir 



force 



17 



19 



20 



21 
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231 ] (i.e., the van der Waals force at larger separations where 



the relativistic retardation becomes essential) and to atom-wall interaction in connection 



with Bose-Einstein condensation 



24 



25|. What is more, in Refs. |26|, |27| the Lifshitz 



theory was extended for the case of an atom (molecule) interacting with a plane surface of a 



2 



uniaxial crystal or with a multi-wall carbon nanotube. The generalization to nanotubes was 



achieved by using the proximity force theorem [28( and the idealization of graphite dielectric 
permittivities which is applicable to multi-wall nanotubes with sufficiently large number of 



walls. In Ref. 



291 ] the Lifshitz-type formulas were obtained for the van der Waals interaction 



between a single layer of graphite (hereafter, graphene) and a material plate, graphene and 
an atom or a molecule, and between a single-wall carbon nanotube and a plate. To achieve 



this goal, graphene was considered in terms of a two-dimensional free electron gas [30]. 
The reflection properties of electromagnetic oscillations on graphene were described by the 
specific boundary conditions imposed on the infinitely thin positively charged plasma sheet, 



31 



291 ] and obtain 



carrying a continuous fluid with some mass and negative charge density 

In the present paper we use the same model of graphene as in Ref. 
Lifshitz-type formulas for the van der Waals interaction between a microparticle (an atom 
or a molecule) and a single- wall carbon nanotube. Both the free energy and interaction force 
at arbitrary temperature are considered. As an example, we have calculated the van der 
Waals free energy and force as functions of separation in the configuration of a hydrogen 
atom and a molecule in close proximity to a single-wall carbon nanotube. Specifically, the 
values of the van der Waals coefficient are determined at different separations. Comparison 
studies of the van der Waals interaction of a microparticle with multi- and single-wall carbon 
nanotubes are performed. If we consider the van der Waals interaction of an atom or a 
molecule with a nanotube having only two or three walls, the idealization of the dielectric 
permittivity of graphite is demonstrated to be already applicable if separation distance 
between a microparticle and a nanotube is sufficiently large. 

The paper is organized as follows. In Sec. II we derive the Lifshitz-type formulas for 
the van der Waals free energy and force acting between a microparticle and a single-wall 
carbon nanotube. Sec. Ill is devoted to the numerical computations of the van der Waals 
interaction between a hydrogen atom or a molecule and a single-wall nanotube. In Sec. IV 
the comparison between the cases of single- and multi-wall nanotubes is performed. Sec. V 
contains our conclusions and discussion. 
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II. LIFSHITZ-TYPE FORMULAS FOR THE INTERACTION BETWEEN 
MICRO PARTICLE AND SINGLE-WALL CARBON NANOTUBE 



We begin with the van der Waals interaction of a graphene occupying the xy-plane, 
z = 0, or of a graphite plate of thickness d, with a semispace made of isotropic material. 
The separation distance between the graphene or the boundary plane of a graphite plate 
(labeled by the upper index 1) and the boundary plane of a semispace (labeled by 2) is a. As 



was shown in Refs. 



26 



27 



291 ]. in both cases the free energy of the van der Waals interaction 



per unit area at temperature T in thermal equilibrium is given by the Lifshitz-type formula 



with the properly defined reflection coefficients r, 



(1,2) . 
TM,TE- 
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1, 2, 3, . . . are the Matsub- 



Here, kg is the Boltzmann constant, £j = 2nkBTl/h with / 
ara frequencies along the imaginary frequency axis (ou = k±_ is the magnitude of the 
wave vector component perpendicular to the z-axis, transverse electric (TE) and transverse 
magnetic (TM) modes are the two independent polarizations of electromagnetic field, and 
qi = (kl + tf/c 2 ) 1/2 . 

From Eq. ([1]), the van der Waals force acting between a graphene or a graphite plate of 
thickness d and a semispace made of isotropic material is: 



F(a,T) 



dF(a, T) 



da 



k B T ^ 
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Now we specify the reflection coefficients. For a semispace made of isotropic material 
(labeled by the upper index 2) they are commonly known 



,(2), 



„(2) 
TM 



£{i&)qi + h 
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where 



kl = dkl + e(z^, (4) 



and s{uj) is the dielectric permittivity of a semispace material. 

If the first body is a graphite plate of thickness d, the reflection coefficients were found 
in Ref. [2d]: 

r (D ( , k , ) _ ft h,\ — eM>Mi)Qi - Ki 



£x{i£i)e z {i£i)q? + k 2 zl + 2qik zt y / 'e x {iii)e z (i&) coth(A^d) 

^te, , ±) = rgUl , , ±) = t , +9? ^;-5 oth( ^ ) . (5) 

Here, s x (ou) = e y {uj) and e^(c<j) are the dielectric permittivities of graphite in the x, y and z 
directions, respectively, and 



Ki = \lk 2 ± + e x (iCi)^, k zl = ^jkl + eMi)^- (6) 



If the first body is graphene, the idealization of the frequency-dependent dielectric per 



mittivity cannot be used. In this case the reflection coefficients can be found [29|, |3l|, |32|] by 



modeling graphene as a two-dimensional plasma sheet carrying a negatively charged fluid of 
7r-electrons. For the hexagonal structure of carbon layers there is one 7r-electron per atom 



331 ] resulting in two 7r-electrons per one hexagonal cell. This leads to the following values 



for the density of 7r-electrons and the wave number of the sheet: 

4 TiC^ 

n = — —, K = 2n — - = 6.75 x 10 5 m" 1 , (7) 

where / = 1.421 A is the side length of a hexagon, e and m are the electron charge and mass, 
respectively. 

Now we are solving the Maxwell equations with the following matching conditions on the 
tangential and normal components of the fields: 

c 2 

E t>2 - E tA = 0, E Z)2 - E Z)l = 2K—V t ■ E t , (8) 

B z , 2 - B z ^ = 0, B %2 - B t ,i = -2iK-j x E t . 

00 

Here j = (0, 0, 1) is the unit vector pointing in z-direction, and all fields depend on time 
through a common factor exp(— iuot). As a result, we arrive to one- dimensional scattering 
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problem in the ^-direction 18) . The solution of this problem leads to the following reflection 



coefficients on a graphene sheet [2f 
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22| : 



TM 



r, 



(i) 

TM, 9 



P (l) 
TE 



(6> k 



± 



r, 



(i) 



,(6,* 



c 2 qiK 
K 



TE, g KV>»±) R+qi 

with the reflection coefficients 



(9) 



Now we use Eqs. (pP) and (j2j) with the reflection coefficients ([3]) and (jU]) to obtain the 
free energy of van der Waals interaction and force acting between an atom (molecule) and 
a single-wall carbon nanotube. For this purpose let us consider an infinite space filled with 
an isotropic substance having a dielectric permittivity e(u) containing an empty cylindrical 
cavity of radius R + a. We introduce the single-wall carbon nanotube of radius R inside this 
cavity (see Fig. 1). In so doing the nanotube is modeled by a cylindrical graphene sheet. 
Then there is a gap of thickness a between the nanotube and the boundary surface of the 
cylindrical cavity. Each element of the nanotube experiences an attractive van der Waals 
force on the source side of this boundary surface. By using the proximity force approximation 
(see Ref. 34J in the case of ideal metals), the free energy of the van der Waals interaction 



between these two co-axial cylinders can be approximately represented in the form 



F c > c (a,T) = 27rL^/R(R + a)J r (a,T). 



(10) 



Here T{a, T) is the free energy per unit area in the configuration of a graphene interacting 
with a plane boundary of an isotropic semispace determined in Eq. flTJ, and L is the length 
of the nanotube, which is supposed to be much larger than R. 

The accuracy of Eq. ffTUl) is rather high. Recently the accuracy of the proximity force 
approximation was investigated on the basis of first principles by comparing the approximate 
results with the exact ones for the configurations of a cylinder or a sphere above plate. It 

that corrections to the results obtained using the proximity 



was shown 



35 



36 



37 



force approximation are less than 0.5a/ R where R is a cylinder or sphere radius. For the 
configuration of two co-axial cylinders the accuracy of the proximity force approximation 
is even much higher. For example, within the separation region < a < R/2, the results 
calculated using Eq. ffTOl) coincide with the exact ones up to 1% for cylinders made of ideal 



metal 



34 



39| (for real materials the accuracy may be different for only a fraction of a 



percent). 
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To come to the case of an atom (or molecule) near a nanotube, we suppose that the 
isotropic substance filling the infinite space is rarefied with some small number N of atoms 
(molecules) per unit volume. We expand the quantity jF c ' c (a, T) on the left-hand side of 
Eq. ffTO]) in powers of N and using the additivity of the first-order term arrive at: 



POO 

^ c ' c (a, T) = N F a > c (z, T)2ir(R + z)Ldz + 0(N 2 ). 

J a 



(11) 



Here, J-" a,c (z, T) is (yet not found) the free energy of the van der Waals interaction of a single 
atom (molecule) belonging to an isotropic substance with a single-wall carbon nanotube (z 
is measured from the external surface of nanotube in the direction perpendicular to it). 
Differentiation of both (negative) sides of Eq. ( TTTT) with respect to a leads to: 



ar c - c (q,r) 

da 



2n(R + a)LNJ ra > c (a, T) + 0(iV 2 ). 



(12) 



This negative derivative can be found also by the differentiation of both sides of Eq. ( flUf) : 

&F c > c (a,T) 



da 

x 



2nL^R(R + 



(13) 



-T{a,T) + F(a,T) 



2(R + a) 

where the force F(a, T) acting between a graphene and a semispace was defined in Eq. 
with the reflection coefficients O) and (EI). 



We expand the dielectric permittivity of a rarefied substance in powers of N 40j 

e(i&) = l + 4vra(^)iV + 0(iV 2 ), 



(14) 



where a(u) is the dynamic polarizability of an atom or a molecule of this substance. Using 
Eq. flH]), we expand also the reflection coefficients rvj^ TE (£z, k±) in Eq. (j3J) in powers of N 
and substitute the obtained expressions in Eqs. and (J2]). Thereafter, with the help of 
Eqs. (1131) and ffl2l in the limit iV — > 0, the desired expression for the van der Waals free 
energy in the configuration of a microparticle interacting with a single-wall carbon nanotube 
is obtained (see Ref. [26( devoted to multi-wall nanotubes for the details of perturbation 
expansion and limiting procedure): 



F a > c (a,T) = -k B T 



R 



R + a 



1=0 



x 



k i dk i e- 2aqi 



A(R + a)_ 



(15) 



g 2 c 2 



In the limiting case R — > oo, Eq. (fl5j) coincides with a known result for the free energy of a 



microparticle near a plane surface of graphene 



291 ] . If, instead of graphene, the cylindrical 



graphite shell of thickness d is considered, Eq. (I15p describes 26] the interaction of a mi- 



croparticle with such a shell if the reflection coefficients r, 



(i) 

TM,TE, 9 



(6) fcjj are replaced with 



P (i) 

TM,TE,d 



(6; k±) defined in Eq. (J5J). 



The van der Waals force acting between a microparticle and a single-wall carbon nanotube 
is obtained as a minus derivative of Eq. ffl5l) with respect to separation distance 



x 



1? 



2«f- 



8(i? + a) 



(16) 



g 2 c 2 



TE, 9 



r (l) 

TM, 9 



By replacing the reflection coefficients ^tm,te, 9 (^' ^-0 with TEd (£j, we return to the 
ex pre ssion for the force between a microparticle and a graphite cylindrical shell of thickness 

d y. 



III. CALCULATION OF THE VAN DER WAALS INTERACTION BETWEEN 
HYDROGEN ATOM OR MOLECULE AND SINGLE- WALL CARBON 
NANOTUBE 

For the purpose of numerical computations it is convenient to introduce the nondimen- 
sional variables under the integrals 

V = 2aqi, Ci = = — (17) 

C UJ C 

and present the free energy (Tl5|) and force (TT6|) in the form 

^ c (a,T) = -^%^, F^(a,T) = -^%^, (18) 
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where 
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8 \l R+a\2(R+a) 



a(0) 



x poo 

1=1 J & 



C F (a,T) 



2{R + a) 

£2 , 

3(2 J R + 3a)(2i? + a 



(19) 



(Ci,y) 



R 
R + a 

OO 



v poo 

+ ^2 a ( i &) / dyye~ y 
i=i •'Ci 



2(,R + a) 2 
3a 2 



a(0) 



x 



4(i? + a) 2 
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(20) 



(i) 

TM, 9 



In terms of new variables (11711 the reflection coefficients fl9l) take the form 



r (i) 

TM,j 



r TE, (0,2/) = r^lJy) 



(Ci,y) 



2yaK 



2yaK + Cr 



(21) 



(i) 



2aK 



TE,g\St) y,/ TE,<?\»7 2rf + y 

To perform computations using Eqs. (fl9j) and (|20|) one needs some expressions for the 



26 



4l|, 



atomic and molecular dynamic polarizabilities of hydrogen. As was shown in Refs 
for the calculation of van der Waals interaction the polarizabilities can be represented with 
sufficient precision in the framework of the single-oscillator model, 

9a 



a(i£i) = a m (i£i) 



9m 



for a hydrogen atom and a molecule, respectively Here, g a 



(22) 



a a (0)u; 2 with the static 



atomic polarizability of hydrogen equal to a a (0) = 4.50 a.u. and the characteristic frequency 
u a = 11.65 eV 4^]. For a hydrogen molecule it holds g m = a m (0)u;^ with a m (0) = 5.439 a.u. 
and uj m = 14.09 eV [42j. Note that when substituting polarizabilities in Eqs. ( Tl9|) . ff20|) . they 
should be expressed in cubic meters (la.u. of polarizability is equal to 1.482 x 10~ 31 m 3 ). 

The computations were performed for single- wall carbon nanotubes with radia R = 2, 3, 4 
and 5 nm at T = 300 K. As an example, in Fig. 2 we plot the van der Waals coefficient (fl9l) as 
a function of separation in the region from 1 to 3 nm for the nanotube with R = 5 nm. Solid 
line 1 labels a nanotube interaction with a hydrogen atom and solid line 2 with a hydrogen 



molecule. For comparison, in the same figure the previously computed [29J van der Waals 
coefficients for the interaction of a hydrogen atom (the dashed line 1) and molecule (the 
dashed line 2) with a plane graphene sheet are included. As is seen in Fig. 2 (the solid lines 
1 and 2), at all separations the van der Waals coefficient of a molecule-nanotube interaction 
is larger than of an atom-nanotube interaction. At the same time, the van der Waals 
coefficients for the interaction of a hydrogen atom or a molecule with graphene (the dashed 
lines 1 and 2, respectively) are larger than the respective coefficients for the interaction with 
a single- wall carbon nanotube. This excess is of about 30%. Note that the obtained results 
practically do not depend on temperature in the temperature region from to 300 K. 

In Table I we present a few computational results for the van der Waals coefficient C3 
in the case of hydrogen atom (columns 2-5) and molecule (columns 6-9) interacting with 
single-wall carbon nanotubes of different radia. From Table I it is seen that the van der 
Waals coefficient C3 is monotonously increasing with the increase of nanotube radius. As in 
Fig. 2, the separation distance between an atom (molecule) and a nanotube varies from 1 
to 3 nm (recall that at shorter separations some additional forces of chemical and exchange 
nature should be taken into account; at larger separations the accuracy of the obtained 
theoretical expressions may decrease). 

In Fig. 3 we plot the coefficient CV for the van der Waals force defined in Eq. (1201) acting 
between a hydrogen atom (the solid line 1) or a molecule (the solid line 2) interacting with 
the single-wall carbon nanotube of R = 5nm radius. As it holds for the coefficient C3, the 
values of the coefficient C f for a molecule are larger than for an atom at all separations under 
consideration. Dotted lines in Fig. 3 (labeled 1 for an atom and 2 for a molecule) represent 
the results obtained on the assumption that Cp = 3C3, i.e., on the assumption that the van 
der Waals coefficient C3 =const and does not depend on separation distance. As is seen 
in Fig. 3, the differences between solid and dotted lines are of about 15-20%. Thus, the 
dependence of the van der Waals coefficients on separation in atom (molecule) — nanotube 
interaction is essential for obtaining the computational results of high precision. 

In Table II the computational results for the van der Waals force coefficient Cp are 
presented. Columns 2-5 are related to the case of hydrogen atom and columns 6-9 to 
hydrogen molecule interacting with single-wall carbon nanotubes of radia R = 2, 3, 4 and 
5 nm. As in Table I, the separation distance varies from 1 to 3 nm. From Table II it is seen 
that the magnitudes of the coefficient Cp increase with the increase of a nanotube radius 
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and decrease with the increase of separation distance. The respective magnitudes of Cf for 
a hydrogen molecule are larger than for an atom for nanotubes of different radia. 



IV. COMPARISON OF INTERACTIONS BETWEEN HYDROGEN ATOM 
(MOLECULE) AND CARBON NANOTUBE IN THE CASES OF SINGLE- AND 
MULTI-WALL NANOTUBES 

In this section we compare the van der Waals coefficients for the interaction of hydrogen 
atom or molecule with single-wall and multi-wall carbon nanotubes. This permits us to de- 
termine how thick should be a multi-wall nanotube in order the idealization of the dielectric 
permittivities of graphite be applicable. We will also quantitatively compare the cases of 
multi-wall and single-wall carbon nanotubes depending on the separation distance between 
a nanotube and a microparticle. 

The van der Waals coefficient of the interaction between a single-wall nanotube and 
a microparticle is given by Eq. ffl9l) with reflection coefficients (j21i . To obtain the van 
der Waals coefficient of the interaction with a multi-wall nanotube, one should replace the 
reflection coefficients ?tm te g m Eq- ( fl9l) with the coefficients v xm,te, cz defined in Eq. (J5]) 
j^ . In terms of dimensionless variables (fTTl) the latter take the form 

e x ie zl y 2 - f 2 z (y,Ci) 



e x ie z iy 2 + f 2 (y, Q) + 2^e xi e z iyf z {y, Q) coth [f z (y, G)d/(2a)] 



r w (c, v ) = fx(y,0)-y 2 (23) 

TE^V) y2 + f 2 {yXl)+2 y fx{yXl)coth[fx{yXl)d/{2a) y K ) 

where the following notations are introduced 

e«ri = e x (iCiu c ), e z i = e z (iQu c ), (24) 

f z (y, 0) = y 2 + (?{e z i - l), f x (y, 0) = y 2 + CK^i - i)- 

The computational results for the van der Waals coefficient C3 of atom-nanotube interac- 
tion are presented in Fig. 4 by the solid dots 1, 2, 3 (the single-wall carbon nanotube of 5 nm 
radius at separations a = 1, 2 and 3nm to a hydrogen atom, respectively) and by the solid 
dots connected with solid lines for the multi-wall carbon nanotubes of 5 nm external radius 
with n — 1, 2, 3, 4 and 5 walls, respectively. For a single-wall nanotube computations were 
performed using Eqs. ( |T9|) . ( 1211 and for a multi-wall nanotubes using Eqs. ( 1T9|) . ( 1231) . In 
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both cases the dynamic atomic polarizability a a is taken in Eq. fl22|) and for a multi-wall 
nanotube d = 3A(n — 1) A. As is seen in Fig. 4, the van der Waals coefficient computed by the 
formula for a multi-wall nanotube with only one wall does not coincide with the coefficient 
computed by the formula for a single-wall nanotube. This is expected result because the 
reflection coefficients in Eq. (I23p approach zero when the nanotube thickness d vanishes. At 
the same time, as Fig. 4 suggests, the van der Waals coefficients for a multi-wall nanotube 
with n = 3 walls at 1 nm from an atom and for a multi-wall nanotube with n = 2 walls at 
a separation 2 or 3nm from an atom are in the proper proportion to the coefficients com- 



puted for a single-wa 
was obtained in Ref. 



1 nanotube using the reflection coefficients ( l2Ti) (the analogical result 



29| for the interaction of a graphene and a graphite plate with an Au 



semispace). This permits us to conclude that the macroscopic concept of graphite dielectric 



permittivity used in Ref. 26( to describe multi-wall carbon nanotubes is already applicable 
for nanotubes containing only two or three walls depending on separation distance between 
a nanotube and an atom. 

Fig. 5 contains the same information, as in Fig. 4, but for a hydrogen molecule. It is seen 
that the approximation of graphite dielectric permittivities is also well applicable for multi- 
wall nanotube with 3 walls spaced 1 nm apart from a molecule and for two-wall nanotube 
at a separation of 2 or 3 nm from a molecule. The values of the van der Waals coefficient for 
a molecule are, however, larger than for an atom in accordance to Sec. III. The single-wall 
nanotube cannot be described using the macroscopic concept of dielectric permittivity but 
as a cylindrical plasma sheet with appropriate boundary conditions. 

V. CONCLUSIONS AND DISCUSSION 

In the foregoing, we have obtained the Lifshitz-type formulas describing the free energy 
and force of the van der Waals interaction of a hydrogen atom or a molecule with a single- 
wall carbon nanotube. The nanotube is represented using the description of a cylindrical 
graphene sheet in terms of two-dimensional plasma carrying free 7r-electrons with appropri- 
ate boundary conditions imposed on the electromagnetic field. Recently this approach was 
used to describe the interaction of a flat graphene sheet with a material plate or a micropar- 



ticle, and a material plate with a single- wall carbon nanotube 29]. The formulas obtained 



here were applied to calculate the van der Waals free energy and force in the configuration 
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of a hydrogen atom or a molecule interacting with single- wall carbon nanotubes of different 
radia. The comparison with the previously developed Lifshitz-type formalism applicable to 
multi-wall carbon nanotubes 26| have allowed to determine the application limits of the 
macroscopic concept of graphite dielectric permittivity. The latter was shown to be appli- 
cable even to two- or three-wall nanotubes depending on the separation distance between a 
nanotube and an atom (molecule). 

The developed formalism suggests some advantages in comparison to phenomenological 
and semi-phenomenological approaches used to calculate the van der Waals interaction. As 



was recently noted in Ref. [Ill ], the phenomenological method of pairwise summation of 
the intermolecular van der Waals potentials can lead to even qualitatively wrong results for 
the nonretarded van der Waals interaction between metallic nanotubes and nanowires (in 
Ref. [43I it was emphasized that the pairwise summation of intermolecular potentials leads 
to correct dependence of the van der Waals potential and force on separation distance only 
for bulk three-dimensional bodies). 

The possibility to use the Lifshitz theory of dispersion forces for the investigation of 
the van der Waals interaction between hydrogen atoms (molecules) and single-wall carbon 
nanotubes demonstrated in the paper may be helpful for the problem of hydrogen storage. 
According to recent review on the subject ^4], there are conceptual possibilities to cre- 
ate carbon nanostructures capable of absorbing more than 10 mass % of hydrogen. The 
resolution of this fundamental problem requires a detailed investigation of the microscopic 
mechanisms of the interaction between hydrogen and graphite. To attain these ends the 
above approach should be extended to distances below 1 nm. At so small separations there 
may be attractive chemical forces and short-range repulsive forces of exchange nature. These 
forces depend on atomic structure and cannot be taken into account by means of the bound- 
ary conditions. Some approximate description is possible by combining the energy of van 



der Waals attraction with phenomenological potentials 451 ]. 
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FIG. 1: Schematic of the cylindical graphene sheet of radius R which is concentrically placed into 
a cylindrical cavity of radius R + a in the infinite space filled with an isotropic substance. 
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FIG. 2: The van der Waals coefficient as a function of separation for the interaction of a hydrogen 
atom (lines labeled 1) or a molecule (lines labeled 2) with the single-wall carbon nanotube of 
R = 5nm radius (solid lines) and with a plane graphene sheet (dashed lines). 
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FIG. 3: The coefficient Cpasa function of separation for the interaction of a hydrogen atom (lines 
labeled 1) or a molecule (lines labeled 2) with the single-wall carbon nanotube of R = 5nm radius 
(solid lines). Dotted lines are drown under the assumption that Cp = 3C3, i.e., that the van der 
Waals coefficient C3 does not depend on separation. 
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FIG. 4: The van der Waals coefficient as a function of the number of walls for the interaction of 
a hydrogen atom with the multi-wall carbon nanotube of R = 5nm external radius (solid dots 
connected with solid lines) and with a single-wall carbon nanotube of the same radius (solid dots 
1, 2, 3) spaced at 1, 2 and 3nm from the atom, respectively. 
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FIG. 5: The van der Waals coefficient as a function of the number of walls for the interaction of 
a hydrogen molecule with the multi-wall carbon nanotube of R = 5 nm external radius (solid dots 
connected with solid lines) and with a single-wall carbon nanotube of the same radius (solid dots 
1, 2, 3) spaced at 1, 2 and 3nm from the molecule, respectively. 
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TABLE I: The van der Waals coefficient as a function of separation for the interaction of a hydrogen 
atom or a molecule with single-wall carbon nanotubes of different radia. 



C 3 (a.u.) 



a 

(nm) 




hydrogen atom 






hydrogen 


molecule 




R = 2nm 


R = 3 nm 


R = 4nm 


R = 5 nm 


R = 2nm 


R = 3nm 


R = 4nm 


R = 5 nm 


1.0 


0.0503 


0.0544 


0.0569 


0.0585 


0.0664 


0.0718 


0.0750 


0.0771 


1.2 


0.0460 


0.0502 


0.0528 


0.0545 


0.0604 


0.0659 


0.0692 


0.0715 


1.4 


0.0424 


0.0466 


0.0493 


0.0511 


0.0554 


0.0610 


0.0644 


0.0668 


1.6 


0.0393 


0436 


0463 


0.0482 


0.0513 


0.0568 


0.0603 


0.0627 


1.8 


0.0367 


0.0410 


0.0437 


0.0456 


0.0478 


0532 


0568 


0592 


2.0 


0.0345 


0.0387 


0.0414 


0.0433 


0.0477 


0.0501 


0.0536 


0.0561 


2.2 


0.0325 


0.0366 


0.0394 


0.0413 


0.0420 


0.0474 


0.0509 


0.0534 


2.4 


0.0307 


0.0348 


0.0375 


0.0395 


0.0397 


0.0449 


0.0484 


0.0509 


2.6 


0.0292 


0.0332 


0.0358 


0.0378 


0.0376 


0.0427 


0.0462 


0.0487 


2.8 


0.0277 


0.0317 


0.0343 


0.0363 


0.0357 


0.0407 


0.0442 


0.0467 


3.0 


0.0265 


0.0303 


0.0330 


0.0349 


0.0340 


0.0389 


0.0423 


0.0448 



24 



TABLE II: The coefficient Cf as a, function of separation for the van der Waals force acting between 
a hydrogen atom or a molecule and single-wall carbon nanotubes of different radia. 



Cf (a.u.) 



a hydrogen atom hydrogen molecule 



(nm) 


R = 


- 2 nm 


R = 3nm 


R = 4nm 


R = 5 nm 


R = 2nm 


R = 3nm 


R = 4nm 


R = 5 nm 


1.0 


0. 


175 


0.186 


0.193 


0.197 


0.232 


0.248 


0.256 


0.262 


1.2 


0. 


162 


0.174 


0.181 


0.185 


0.214 


0.230 


0.239 


0.245 


1.4 


0. 


150 


0.163 


0.170 


0.175 


0.198 


0.214 


0.224 


0.231 


1.6 


0. 


140 


0.153 


0.161 


0.166 


0.184 


0.201 


0.211 


0.218 


1.8 


0. 


132 


0.145 


0.153 


0.158 


0.173 


0.190 


0.200 


0.207 


2.0 


0. 


124 


0.138 


0.146 


0.151 


0.162 


0.180 


0.190 


0.197 


2.2 


0. 


118 


0.131 


0.139 


0.144 


0.154 


0.170 


0.181 


0.189 


2.4 


0. 


112 


0.125 


0.133 


0.139 


0.146 


0.162 


0.173 


0.181 


2.6 


0. 


108 


0.120 


0.128 


0.133 


0.138 


0.155 


0.166 


0.174 


2.8 


0. 


102 


0.115 


0.123 


0.128 


0.132 


0.148 


0.159 


0.167 


3.0 


0.0975 


0.110 


0.119 


0.124 


0.126 


0.142 


0.153 


0.161 
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